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INTRODUCTION 



The 3+1 formulation of general relativity views spacetime as the time history of the 
geometry associated with a spacelike hypersurface S. In the differential-forms (df) version^ 
of 3+1 relativity, the (co)triad 1-form and the "extrinsic curvature 1-form" are taken 
as the basic geometrodynamical variables of interest. [|l[ In this description initial data sets 
(e*^, K^) are constrained by seven point-wise relations, DF versions of the familiar Hamil- 
tonian, momentum, and rotation constraints; and the Arnowitt-Deser-Misner (adm) 
equations of motion for and govern the dynamic evolution of the system. These evolu- 
tion equations preserve the constraints, and, moreover, maintain complete spatial covariance 
by means of an arbitrary shift vector (3^ and an arbitrary skew-symmetric rotation matrix 
0ab (covariance under spatial coordinate transformations is not an issue in the DF language). 
However, as is well-known these equations do not manifest mathematically the propagation 
of physical degrees of freedom along the lightcone; they do not give rise to a physical wave 
equation. Therefore, in 3+1 form relativity (as in its pure-index metric counterpart) gauge 
and physical variables remain intertwined. 

Recently, Choquet-Bruhat and York 0] and Abrahams, Anderson, Choquet-Bruhat and 
York [|] (hereafter, we refer to both citations jointly as aacy) have cast the 3+1 evolution 
equations for metric general relativity in gauge-covariant and causal "first-order symmetric 
hyperbolic form" (fosh form), thereby cleanly separating physical from gauge degrees of 
freedom in the Cauchy problem for general relativity. A key ingredient in their construc- 
tion is a certain wave equation which governs the light-speed propagation of the extrinsic 
curvature tensor Kij. Along a similar line, we construct a related wave equation which, as 
the key equation in a system, describes general relativity. (Here we only consider vacuum, 
although we expect the inclusion of matter to be straightforward enough.) Whereas the 
approach of aacy is based on tensor methods, the present formulation is written solely in 
the language of differential forms. Our approach starts with Sparling's tetrad-dependent 
differential forms, and our wave equation governs the propagation of Sparling's 2-form, 
which in the "time-gauge" is built linearly from the extrinsic curvature 1-form. The tensor- 
index version of our wave equation describes the propagation of (what is essentially) the 
Arnowitt-Deser-Misner (adm) gravitational momentum. As with aacy's "Box K" equa- 
tion, to be proper our DF wave equation requires the implementation of a harmonic time 
slicing (or suitable generalization thereof) in order to handle the time re-parameterization 
invariance of Einstein's theory. In this paper we only derive the described wave equation. 
We do not discuss the possibility of constructing FOSH systems for Einstein's theory related 
to the one given by aacy. (We hope to return to this issue later.) Our goal is to understand 
the geometry of the aacy formalism in the DF language. Although there is no doubt in 
our mind that the original pure-index formulation of AACY is the superior wave-equation 
version of Einstein's theory, one might expect a DF construction to show distinct advantages 
in some theoretical contexts. 



^ Quite a number of reviews of "form relativity" exist. An excellent choice is the discussion found 
in Ref. by Isenberg and Nester. 
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AACY obtain their wave equation for Kij by (i) taking a further "time derivative" of the 
equation of motion for K^j, and then (ii) subtracting from the resulting expression certain 
spatial covariant derivatives of the momentum constraint (which has (3^ as its corresponding 
Lagrange multiplier). Now, the standard DF version of 3+1 relativity is, in essence, a triad 
theory, requiring the existence of three auxiliary phase-space constraints (packaged into 
one name, the "rotation constraint," with corresponding Lagrange multiplier 0ac)- In our 
wave-equation construction, we proceed in a similar fashion to AACY. We also take the 
time derivative of an equation of motion and then subtract from the resulting expression 
certain spatial covariant derivatives of the momentum constraint. However, we find the 
need to also subtract from this time derivative certain spatial covariant derivatives of the 
rotation constraint. To us this seems to be natural generalization of the AACY method, 
necessary to handle the extra gauge freedom present in the triad theory. It is well-known 
that geometrically the rotation constraint is related to a certain piece of the spacetime 
torsion; and, therefore, in this paper we keep track of this very piece. Of course, this has 
consequences for the usual symmetries enjoyed by the Riemann, Ricci, and Einstein tensors 
(and DF versions thereof). However, we wish to stress here that we do not consider Einstein- 
Cartan theory in this paper. "At the end of the day," we set all torsion equal to zero. We 
keep track of torsion in the beginning only as a rigorous way to handle the issue of the 
rotation constraint, which is always present in a triad version of Einstein's theory. 

The organization of this paper is as follows. In a preliminary section, § I, we fix our 
conventions and collect the geometric tools necessary for our central discussion. In particular, 
we quickly review the 3+1 decomposition of spacetime differential forms and introduce a 
certain evolution operator Dq, which may be regarded as a (triad) generalization of the 
evolution operator central to the work of AACY. The reader may wish to simply peruse 
the preliminary section, and if necessary return to it later for more detailed study. § II 
deals with geometry of vacuum spacetimes as described by differential forms, and it lays 
the groundwork for our wave-equation formulation of general relativity. § HI introduces 
Sparling's forms, using them to enact a DF version of the standard 3+1 decomposition of 
spacetime geometry, while § IV combines the results of § II and HI to write down a system 
of DF equations for general relativity, one equation in the system being the promised wave 
equation. We demonstrate the equivalence between our DF system and the original Einstein 
equations in § V. The concluding section provides some "instructions" for rewriting our 
results in tensor-index form. Finally, we have included as appendices a handful of technical 
derivations necessary for our central discussion. 

I. PRELIMINARIES 

A. Notation and Conventions 

Consider a spacetime M., equipped with Lorentz-signature metric g^y, which is foliated 
by a collection of spacelike 3-dimensional slices S. We use S both to denote the foliation and 
to represent a generic slice. Since we are concerned only with local field equations in this 
work, we do not further specify the global properties of S. The timelike, future-pointing, 



3 



unit, hypersurface normal of the S foliation is u'^. Respectively, hij and Kij denote the 
intrinsic metric and extrinsic curvature tensor associated with the S foliation. 

We shall use two types of spacetime frames in this work. One is the foliation-adapted 
"quasi-coordinate" frame (with coframe e'^) used by aacy 

6° = dt eo = d/dt - (3^d/dx^ 

(1.1) 

= dx^ + (3^dt ej = d/dx^ . 

Greek letters {//, i^, X, ■ ■ ■} run over 0, 1, 2, 3, but is not t. Here N and are respectively 
the standard lapse function and shift vector associated with the E foliation. Lowercase Latin 
letters {i,j, k, ■ ■ ■}, taking the values 1, 2, 3, represent the space legs of the quasi-coordinate 
frame and are also S coordinate indices. We use do to represent eo as a Pfaff derivative, § 
so for some scalar function / we write dof to mean / — f]^df /dxK 

The second type of frame we consider is the foliation-adapted (or time-gauge) tetrad ca 
(with coframe e^): 

= A^e" e± = N-^cq 

(1.2) 

e c jC"^ e^ ^a'^^j • 

Notice that ej_'^ is the S normal m^. Capital Latin letters {A, B,C,- ■ ■} represent tetrad 
indices and run over ±, 1, 2, 3, while block lowercase letters, typically from the first part of the 
alphabet, {a, b, c, ■ ■ ■} represent S triad indices and run over 1, 2, 3. The set of ea-^ comprise 
an orthonormal triad of vector fields on E. Orthonormality implies that hijCa^eh-^ = S^h, 
with (5ab = = diag{l, 1, 1). We can represent the S metric as hij = S^hG^ j- 



B. 3-|-l Splitting of Forms and the First Cartan Structure Equation 

We may split a general spacetime differential form \1' into pieces normal and tangential 
to the S foliation, 

^ = ii^ + e^A^x, (1.3) 

where the i± = iu is the "inner derivative" with u, also known as the vector-form "hook" 
with u. The hook of a spatial form with the timelike S normal u = ej_ vanishes. Both 
and \E'x are spatial forms. Let us take a special look at spacetime forms which happen to be 
spatial to start with. For the sake of definiteness only, consider a spatial 2-form ip = ^kp. In 
the index notation and with the quasi-coordinate basis, the only non- vanishing components 
of '^p^ll' are ipij, so that ip = ^ ■ Note that, while ip is purely spatial, it is still a 

spacetime form; we have not pulled it back to a particular E slice. This is an advantage 
gained by using the quasi-coordinate frame: all forms (in fact, all tensors) which are parallel 
to the E foliation in spacetime may be represented with E indices in the index notation. As 
another example, take = "e'^, which has non- vanishing components e'^ j. 



4 



Now suppose \l/ is a 3-form or less. We have the following decomposition of d\l/ with 
respect to the S foliation: 



= c/ii^ + A^-^e^ A [So'i^ - d{N^±)\ , (1.4) 

where d is the Ai exterior derivative and d is the S exterior derivative. Also, our "time- 
evolution operator" is do = £eo, with £ representing Lie differentiation.^ Note that do 
is not the Pfaff derivative do = d/dt — jS^d/dx^, although it agrees with do when acting 
on scalars. In the index notation, the action of (9o on a S form, for example ip\ = ipiC^x, 
is given by £eo'^x = ^^xi'^Pi — Lptpi), with L denoting intrinsic S Lie differentiation. Thus, 
while working (in the index notation) with E tensors carrying spatial indices k,- ■ ■), it 
is consistent to set do = d/dt — Lp. 

Applying the exterior derivative decomposition ( p..4|) twice to the formula (^ip = 0, where 
now ip = ^kp is an arbitrary spatial form, one proves that the following commutator vanishes: 

[do, d]4' = dodip — ddoip = . (1.5) 

In index notation (|1.5|) can be translated into the statement that do commutes with dk 
partial differentiation. However, for the Pfaff derivative [c^o, dk] 7^ 0. 
Consider the first structure equation of Cartan, 

= = de^ + B A , (1.6) 

where = ^e^ xT^ fj,ue^ A e'^ is the torsion 2-form, the T^b are the 1-forms specifying the 
spacetime connection with respect to the time-gauge tetrad (|1.2|), and V is the spacetime 
exterior covariant derivative. We shall assume that T-^ij is the only piece of the spacetime 
torsion tensor which is not manifestly zero; hence we only need to consider the T = "T-*- 
piece of the torsion 2-form. 

The II projections of the (±, a) values of (|1.6|) tell us that 



= De"" = de^ + uj\ A (1.7a) 

T = K^^e^ A , ([T^b) 

where = 'T'^b are the S triad connection 1-forms and D is the S exterior covariant 
derivative. The top equation says that the intrinsic connection on each slice S is torsion- 
free. The second equation follows because (i) e-*- A de^ = and (ii) the triad components of 
the S extrinsic curvature tensor are defined in the time-gauge by K^h = —^a±h- Therefore, 
setting T equal to zero enforces the symmetry of Kij. Taking the i± inner derivative of ( p.. 61) , 
one finds that 



^It is precisely by choosing to take the Lie derivative along cq'^, that we may use spacetime £ in 
( |1.4| ) in place of the "intrinsic E Lie derivative generaUzed to take care of derivatives in directions 
off of the S hyper surf aces" discussed in Ref. Q. Had we chosen to take the Lie derivative along 
u^, as is often done, then we would have needed to introduce such a generahzed E Lie derivative. 
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(1.8a) 



= -ab + r^b±, (&) 

with ttb = eb[logA^]. It follows that the are the triad components of the acceleration of 
the E foliation normal u^. |[l|] Upon inspection of ( |1.8| a), we define a generalized version Dq 
of the evolution operator do by 

^oe^ = 9oe" + r%oe\ (1.9) 
and will shortly extend its action to other forms and tensors. 



C. Triad-covariant Evolution Operator 



The advantage gained by using the simple evolution operator Oq is well-known. |T2 
Differentiation by do preserves the spatial character of any spatial form or spatial tensor. 
Moreover, do defines a good "time axis" as it is orthogonal to the spacelike S foliation. Our 
generalization of do, the degree-zero evolution operator Dq introduced in ( |1.9| ), is covariant 
with respect to triad or "internal" rotations. Now, in fact, the canonical formulation of triad 
gravity identifies the anti-commuting Lagrange parameter (j)'^ b associated with the rotation 
constraint as — F'^bo = — -?Vr^b±- The connection coefficients F^ao describe the rotation of 
the triad as it is parallel transported along the integral curves of cq'^ = Nu^. We generalize 
the action of Dq to other triad tensor- valued forms in the obvious way. For example, consider 
a triad tensor-valued 0-form fj". On such an object, the direct way of expressing the action 
of Do would be 

DofJ" = DofJ" ^ dofj" - /c"F\o + /aT^o • (1.10) 

Notice that on 0-forms Do reduces to Dq. If we had taken fj° as a tensor- valued 1-form 
fa°kG^) then we would have had to use the operator do in the first term on the right-most 
side of ( |1.10|) . In other words, the "hat" does not see triad indices, rather it sees whether or 



not the differentiated object is a form of degree higher than zero. Equation ( |1.8| a) ensures 
that we have 

Doe'' = -NK'', (1.11) 

with the extrinsic curvature appearing as a one- form K'^ = K'^jc^ , as the evolution rule for 
the (co)triad. 

Let us address some technical points concerning how to use the operator Do in index 
notation. First, when working in index notation, assume the following evolution rules for 
the cotriad and triad: 

Doe", = doc^j + e^,F\o = -NK^j (1.12a) 
Doe,^ = doej - eb^F\o = iVeb^ir\ • (|FT^ ) 
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Second, assume that Dq obeys the Leibnitz rule and reduces to do when acting on tensors 
without triad indices. As an example of how to proceed in index notation, consider again 
our triad tensor-valued 0-form, but assume that its triad components fg^ = fi^e^^e^j are 
obtained by "soldering" all of the free indices of a S tensor fi^. One can also obtain the 
action of Dq on f^^ by letting it act on fi^Cg^'^e^j and then using the Leibnitz rule. The result 
is 

Dofe^ = ea'e^^o/,^' + Nf.'^K^ - Nfg,'K\ . (L13) 



Of course, ( |1.10D and ( |1.13|) agree. By first writing 6ah = Kj^a^^h^ and eabc = ^ijk^a^h^^c 



where eijk is the S permutation tensor fixed by 6123 = one can use the result d^hij — 
—2NK(ij^ to show that Dq annihilates both (5ab and eabc- That Dq kills these particular triad 
tensors can also be shown with formulae like (|L1CI| ) along with the fact that the Fabo are 
antisymmetric in their triad indices. 

We have already seen that do and d commute when acting on purely spatial forms. What 
about Dq and Dl {D is the S exterior covariant derivative.) Appendix B proves that for 
any triad-valued S form tp^ = ^kp^, we have the following commutation rule: 

[Dq, D]ij^ = {DqD - DDQ)ij^ = N (3?^ b± + a^K^ - a^K^) A , (L14) 

where 3fJ^b± = 'i±^^h and the SfJ'^b comprise the triad block of the matrix-valued curvature 
two-form of Cartan, = ^^"^ Bfiu^^ A e^. Here the Bfiu = €.'^a€-B'^^°'i3fiu are the mixed 
components of the spacetime Riemann tensor. 



II. EINSTEIN GRAVITY IN TERMS OF DIFFERENTIAL FORMS 

Working within the Cartan framework of tensor-valued differential forms, in this section 
we collect a handful of standard results concerning the geometry of vacuum Einstein space- 
times. We carefully examine the role played by the residual torsion in our formalism, in 
particular deriving several identities which will be used later on to demonstrate that the 
residual torsion may safely removed. 



A. Einstein 3-form 

Consider the Einstein tensor G/ja = e^^G^a written with respect to a mixed basis. Here 
and it what follows we pay strict attention to index ordering. Since we are keeping track of 
certain components of the spacetime torsion, the Einstein tensor is not symmetric a priori 
in our formalism. We will work with the Einstein 3-form 

GA = lG^Ae\fsxe'' Ae^ Ae\ (2.1) 

Clearly, *Ga = G^ag^, where * denotes the Hodge-duality defined by the metric g^i,. In 
the A4 permutation tensor ex^^a is chosen so that 60123 = V—Q- As described in the 
preliminary section, decompose the Einstein 3-form as 
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Ga = "Ga + a Ga± 



(2.2) 



Letting -k denote the Hodge-duality defined by the S metric hij, we have the relations 
*Ga± = —GjAG^ and *"Ga = —G±a- When working with the * duality in index notation, 
the S permutation tensor eijk is the appropriate one to use. 

Let us write down DF equations which capture the antisymmetric pieces of the Einstein 
tensor. These will be related to the residual torsion Tij in the next subsection. First, one 
immediately finds that 

Gije' A = e'^ A *Ga± . (2.3) 

Now turn to the antisymmetry in the cross components G±j — Gj±. We have already seen 
that *"Ga = —G±g,. To get the other cross term look at G^_l, which obeys '^Gj,^ = —Gj±e^. 
Therefore, the 1-form expression 

(G^, - G,^)e' = -(*iiGa)e^ + '^G^x (2.4) 
captures the cross-term antisymmetry. 



B. Bianchi Identity for the Torsion 2-forni 



In this subsection we shall relate the residual torsion T to the antisymmetric pieces of 
the Einstein tensor. To do this, we apply the form decompositions (1.3, L^) to the Bianchi 
identity for the torsion. 



VT^ = dT^ + A = 3? 



(2.5) 



subject to the assumption that both T'^ = and i±T^ = 0. 



The _L value of (|2.5|) is then 



dT = A . 



(2.6) 



The " projection of (|2.6| ) leads to an equation for the E expression dT which is of no concern to 
us here. However, taking the i± inner derivative of ( |2.6| ) and using the form decompositions 
( |1.3| , |1.4|) , we obtain the following identity: 



N-'doT = 3ft^a±be^ A e"" 



(2.7) 



In order to re-express ( p.7|) in terms of the Einstein 3-form, first use the relations^] 3fJ"'"a±b 
ea*eb-'3ftij — ^^s,ch and 3fJ[jj] = G[ij] to rewrite the preceding identity as 



Gi.e' A = -N~^doT + ^^chc'' A 



(2.8) 



In these formulae represents the spatial components of the spacetime Ricci tensor. 
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Next, plug the Gauss-Codazzi-Mainardi (gcm) equation (|Al|a) into ( |2.8| ), use the fact that 
the S Ricci tensor Rij is symmetric (since we have spatial torsion T^jk = 0), and appeal to 
(p.3|). These steps lead to a key identity 

doT = NHT + NK^ A *(e^ A *T) - Ne"" A *Ga± , (g^) 



which gives the "time derivative" d^T of the our residual torsion. This important equation 
establishes that, given the vanishing of the torsion Tj^ pulled back to some initial Cauchy 
surface, no torsion can develop for vacuum spacetimes characterized by a vanishing Einstein 
3-form. 

We now turn to the triad values of the Cartan-Bianchi identity ( p.5| ) for the torsion, 
which for our restricted torsion take the following form: 

AT = 3fJ^BAe^. (2.10) 



Using the GCM equation ([Al|a) for the 3 + 1 splitting of the triad components of the space- 
time Riemann tensor, one finds that the " projection ( |2.10| ) is an identity which is trivially 
satisfied. However, taking the z_l inner derivative of ( p. 10 ), one obtains 

a^T = V A + afJ'^x , (2.11) 
which after a few more manipulations yields 

(*iiG^)e*-Gx± = aA^T. (2.12) 
Here and in what follows we make extensive use of the following the common short-hands: 

= eabce^ (2.13a) 



e 



* _ 1 ^ J) 



.a - leabce^ A . (|27r3| b) 

One should understand that *"G^ is a 0-form. Looking back at (|2.3| ), we see that this identity 
relates the torsion 2-form to the cross-term antisymmetry of the Einstein tensor. In index 
notation, the equation 

G,^-G^i = T,,a^ (2.14) 

captures this antisymmetry. 



III. SPARLING FORMS AND THE 3 + 1 EINSTEIN EQUATIONS 

Sparling's differential forms have played a central role in several recent developments in 
general relativity, [1^ and the tetrad versions 0,^ of these forms also play a principal role 



in our formalism. We begin this section by introducing the tetrad Sparling forms and by re- 
viewing their remarkable properties. In particular, we recall that the Sparling forms provide 
a handy expression for the Einstein 3-form. Meshing quite nicely with our decomposition 
(O) of exact spacetime forms, this handy expression leads to a 3 + 1 decomposition of the 



Einstein 3-form which is particularly useful in our construction of a nonlinear wave equation 
for general relativity. We consider the special case when the (tetrad-dependent) Sparling 
forms are determined by a time-gauge tetrad (|1.2| ). Although the results of this section are 
standard ones, it would appear that we perform the 3+1 decomposition in a novel way. 
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A. Basic Properties and Time-gauge Expressions 



We start with the Sparling relation, which expresses the Einstein 3-form 

= dcTA -TA- \eABCDT'' A T^^ (3.1) 

in terms of the torsion 2-form, the tetrad connection coefficients, and the tetrad-dependent 
Sparhng 2-forms and 3-forms, 

a^^-lr^^Ae^ (3.2) 

= i (r^ ^ A A e*j,cD - A T^^ A e\j,c) ■ (3-3) 

Here e*A^Q = eAscD^^ ■ Notice that, due to the way in which the tetrad connection coef- 
ficients appear in the definitions above, the Sparhng forms do not behave homogeneously 
under tetrad transformations. However, while individually neither dcr^, r^, nor the final 
torsional term in the Sparling relation ( |3.1|) behave homogeneously, the sum of terms on 
the right-hand side of ( |3.1D must behave homogeneously under tetrad transformations, since 
this sum is the Einstein 3-form. 

We shall have need to consider only the Sparling forms {a^, a^, rj_, Ta} determined by a 
time-gauge tetrad (|1.2|) . Consider first the " and projections of the 2-forms {(J_l, cTa}. One 
quickly finds the following list: 

iiaa = el^ A if ^ (3.4a) 
aa± = c^a + a^e^b dUb) 
"cTx = ^a A e'^ (0c) 
= iV-iraoe% dOd) 

where for notational convenience we define tUa = — |eabci^^'^- We make the same definitions 
for Fao and the triad- valued curvature two- form below in ( |3.5| c). Notice that the forms "(Ta 
and cTa± — C(Ja, behave homogeneously under triad transformations, and hence that it makes 
sense to consider their exterior covariant derivatives, -Diba and Dio^^ Evidently, the 

expressions for a^^ and 'bx behave inhomogeneously under rotations of the triad. 

The task of finding the explicit time-gauge expressions for the set {"ra,ra±, "r^,Tj_^} is 
somewhat tedious, and to get the expressions in the list below, one must make extensive use 
of Cabc permutation-symbol gymnastics. We find the following set: 

lira = A if ^ A e^, + ieabc^'^^ A T (3.5a) 

llrx = iif ^ A if ^ A - A A e^^ (H^b) 
Tax = -A^-^r^ aoA^' A 6^, \ - du^ - a N 

+a^a;^ a A + kabcif A if ^ - N'^Y^^T (pic) 



"be "T 2 

r^^ = -a^if ^ A - |eabcif " A cu^^ - N'^V^ do^"" A e^^ , (Ipd) 



where in (pr3|c) one finds the S triad-valued curvature 2-form i?a = duj^ + ^Cabci^^ A uj"" . 
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We shall be dealing exclusively with the projected 2-form "(Ta. Therefore, to avoid undue 
notational clutter, let us define the S form 



7ra = -|ilaa=iK"Ae:b (3.6) 

(the factor of — | has been chosen for later convenience) and collect a few results concerning 
it. First, with the definition ( |2.13| b), the evolution rule ( |1 . 1 1| ) for the cotriad shows that 



Doe: = -2iV7ra. (3.7) 

Moreover, defining^ 



TT 



''^ - e'' An^ = liHS"^ - K''^)e* , (3.8) 



we take the rotation- constraint 3-form to be 

In fact, J^^ = , where the are the triad components of the residual torsion 

2-form. Alternatively, the construction 

T = 2e^AVa (3.10) 

provides a direct way of expressing the torsion 2-form. 

We wish to point out that 2-form vTa is essentially the momentum conjugate to the cotriad 
in the Hamiltonian formulation of triad gravity. Indeed, (switching to index notation) 
the E Hodge dual of vTaij is 



8nh-'/^e,,n''J,,,, (3.11) 



where vr^^.^^^ = {lQn)^^\/h(Kh^^ — K'-^) is the standard ADM gravitational momentum. 
Therefore, were we considering the standard Hamiltonian triad formalism, we would have 
the following ClS db canonical pair: 

{(87r)-Va.,,g^*^} , (3.12) 

where e^^^ = \/he^^^e^ k- The standard rotation constraint which arises in Hamiltonian triad 
gravity (also known as the Gauss constraint in the Ashtekar formulation [0) is nothing but 
our J^^ (stripped of a (Px). The vanishing of the rotation constraint is equivalent to the 
vanishing of our residual torsion 2-form T. 



^Note that this is a special definition, and that vr^*^ is not obtained from vr^ via inner differentiation, 
i. e. TTab ^ «bvra. 
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B. 3+1 Decomposition of the Einstein 3- form 



Using the form decompositions ( |1.3[ ) and ( |1.4| ) and recalling that the only non-zero piece 
of the torsion 2-form is T = 'T-*-, we immediately write the Sparling relation ( p.l|) as the 
following _L and " projections: 



G.I =iV-i 



(3.13) 



Now, in fact, we know that as a whole these pieces must be triad-gauge covariant. Indeed, 
insertion of the explicit expressions ( p.4| ) and ( p.5|) into ( |3.13|) yields the following nice 
formulae: 



"Ga = -2Dtx^ 

"G^ = 2VabVr^" - H7i\ + A i?a 

Ga± = -2iV^iDo7ra + N-^el^ A D(iVa'^) - 



(3.14a) 
dUlb) 

(EH) 



where H = and vr'^a = He*. Notice that ( p.l4| a) and ( |3.14] b), respectively, are 3-form 
versions of the standard momentum and Hamiltonian constraints.^ We could effortlessly 
obtain the 3+1 expression for G^^ with ( |3.14| a) and our previously derived result ( ^^.12] ). 



IV. WAVE EQUATION FOR GENERAL RELATIVITY 

Let us quickly take stock of the situation thus far. As a statement of the standard 3+1 
version of Einstein's theory in vacuo, we take Eq. ( |3.7| ) (in practice, the equation of definition 
for TTa), the vanishing of Eq. ( p.l4| c) (the equation of motion for VTa), and the vanishing of 
Eqs. ( |3.9| ), ( ^.14| a), and ( |3.14| b) (rotation, momentum, and Hamiltonian constraints). After 
making a short digression in § IV.A in order to discuss the Laplacian and wave operators we 
shall use in our construction, we introduce in § IV. B a certain triad-valued 2-form Ta whose 
vanishing will, in effect, replace the evolution equation for vTa. We then perform a 3 + 1 
decomposition on the 2-form Ta. This decomposition yields a non-linear wave equation for 
TTa (assuming a harmonic time-slicing) which serves as the key equation in our system of DF 
equations, written down in § IV. C, which describes vacuum general relativity. 



^In in Isenberg and Nester expressed the gravitational constraints in essentially this way. The 
tt" defined in a footnote on page 60 of that reference is, apart from an overall factor of —4, the TTa 
that we use here. Also, their tt"* differs from our tt*^^ by a factor of —2. 
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A. Digression: Laplacian and Wave Operators on S Triad-valued 2-forms 

Define a triad-covariant Laplacian A whose action on triad- valued S 2-forms ipa, = "■i/'a is 

A^, = (*D'D-D'D*)^^. (4.1) 

Clearly this operator reduces to the standard DF Laplacian when acting on ordinary 2-forms. 
Our chief interest is the action of the triad-covariant Laplacian A on the 2-form vTa. A short 
calculation shows that 

AvTa = ea^e^XV™V'"7r,, - n^.R^, - n^mp,,m)el , (4.2) 

where vTy = ^{Hhij — Kij) and, in accord with the notation of aacy, we use to denote 
the torsion- free covariant derivative operator compatible with hij. Note that vr*-^ differs from 
^ADA/ by the density factor {8n)^^^/h. Also, Rij is the S Ricci tensor and Rijki is the S 
Riemann tensor. Clearly, even upon passage to the tensor-index formalism {*D*D — D*D*) 
differs from the tensor V'^Vfc used in aacy. However, provided that equivalence with the 
Einstein equations has been shown, the extra curvature terms buried in the form Laplacian 
can be eliminated in favor of the variables of interest via a definite prescription. The 
ability to eliminate S curvature terms has proven to be crucial in aacy's construction of a 
(gauge-covariant and causal) FOSH system for Einstein's theory. 
Our wave operator isQ 

6n.,^[-N-\Do)' + A]rr,. (4.3) 

Writing nvTa = eJe^^Aijd^, one finds that 

A,, = [-iV"2(4)2 + v'=Vfc]7r,, 

-|- [other terms in VTjj and ^oVTjj (but no second derivatives of i^ij)] , (4-4) 

which more than hints at the fact that the tensor-index formulation of our construction 
involves a wave equation for Ttij. The tensor wave operator in ( [4.4|) is one used by aacy. 



B. Taking Derivatives of the Field Equations 

The aforementioned 2-form Ta is defined by the following equation:[] 

2Ta = ^oGa± - iV*Z}^iiGa " A D^mQ"^) + iVa"(*iiGb)e: - D(iVTa) (4.5) 

(the unimportant factor of 2 sitting before Ta is chosen only for later convenience). Here Ta is 
a certain carefully designed 1-form which is built purely from torsion (and spatial derivatives 



^We could, of course, work instead with the slightly different operator = —{N ^D^)'^ + A, but 
we shall stick with 



^See the Eq. (|5.4D for the tensor-index definition corresponding to Ta. 
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of torsion) whose particular form [given in the appendix Eq. (|C6|)] need not concern us yet. 
Notice that Ta has been constructed by taking a "time derivative" of the equation of motion 
( p.l4| c), with certain spatial derivatives of the momentum and rotation constraints subtracted 



from the result. Let us briefly comment on the origin of the various terms in the definition 
( [4.5| ). In performing the 3 + 1 decomposition of DoG^x, one generates certain Einstein 3-form 
and torsion terms, in addition to other terms which are built with neither Einstein 3-form 
nor the residual torsion. Now, in fact, the subsequent terms in Ta have been tailored to 
exactly cancel these generated terms. However, we note that our expression for Ta is not 
unique, as one may always add terms to Ta which are homogeneous in the Einstein 3-form 
and torsion. Indeed, the fourth term in ( [4.5| ), linear in the undifferentiated Einstein 3-form, 
may be striken from the definition. We have included this term in order to obtain a more 
pleasing result for Ta's 3+1 decomposition. However, this would seem to amount only to a 
matter of taste. 

To obtain the 3+1 decomposition of 2-form (O), perform 3+1 decomposition on the 



very first term -DoGa± in the definition. With the splitting result (|3.14| c), a straightforward 
expansion gives 

DoGax = -2£)o(iV-i^o7ra) + Do[N-'el^ A D(iVa^)] - DoR. 

+2{doH)7T, + 2HDon^ + ^o(2VbcV*e: - H^el - 4Xavr^) . (4.6) 

The second and third terms on the right-hand side of this identity are the ones which require 
some work. The third term, —DqR^, is handled in Appendix C; therefore, focus attention 
on the second term. Using the evolution rule (|1 . 1 1|) for the cotriad and the commutator 
result ( |1.14| ), we write this term as follows: 

bo[N-^el^ A D{Na^)] = -Na^i^'^'^x + 2a[^K'^l) A e^^ + Cabc^^ A D(A^a'=) 

-N-'^{doN)el^ A D{Na^) + N~^el^ A DDo{Na^) . (4.7) 

At this point, we perform some permutation-symbol gymnastics with the first term on the 
right-hand side in order to reachQ 

bo[N~'e:^AD{Na'')] = 

Na^G±^el - |eabcA^(3?'''± + 2a^K') A a + Na^a^iHS'^'' - K'°^)el 

+eabci^'' A D{Na^) - N^\doN)el^ A D{Na°) + A^~^eab A Dbo{Na°) . (4.8) 

Next, with the index evolution rule (|1.12|) we derive the identity 

Do{Na^) = Doie^'diN) = D^{doN) + N^a^K^^ (4.9) 

(on a scalar like OqN, the action of is e°^di), and then insert it into the last term of 
( [4.8|) , expanding the resulting expression. At the same time, we make a substitution into 



^We note that there is another possible route which one could take at this point. Namely, one 
could directly substitute the appendix result (p^ ) into the first term of (4/7). If one intends to 
follow this route, then one should first strike the term A^a''(*"Gb)ea from the definition (|4.5|) of Ta- 
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the second term in ( |4.8| ) with the appendix result ( |C7| ). Combining the result of these 
endeavors with ( [4.6|) and the appendix expression (|C8|) for -Do-Ra, we readily obtain the 
desired expression for Z)oGa_L. Finally, all extrinsic curvature terms in the expression are 
re-expressed in terms of vr*^^ via (|3.8|) . 

With the splitting result for DoGa±, we express the 3+1 decomposition of the 2-form 
(BTSl) as follows: 



Ta = iVDTTa + f/I + f/J + K , (4.10) 

where explicitly one has the source term 

U: = ^o(VbcV*e: - '^H^el - 2Vba7r^) + {doH)n, - 2Na A *D*n, 

-iV-^e^b A DiN'Ha'') + N'^fDoTc, - ^N-'fel^ A /^(iVa^) , (4.11) 
the torsion term 

Ul = l5^bN^[*diN''T)]a A - i(5abiV'a A *D*{N~'e^ A *T) + ^N-'D{N^a/T) , (4.12) 
and the slicing term 

K = liV-ie^b A D(Z}V) • (4.13) 

In these expression, we have defined / = OqN + N^H. Notice that all of the second- 
derivative terms in (|4.1CI| ) involving vTa (other than those found in nvTa) have been collected 



together in the slicing term ( |4.13| ). Also, make note of the torsion appearing in U^. Here we 



consider T as a short-hand for 26*^ A Va [cf. Eq. ( p.lO| )]; and, therefore, in effect we have the 
combined expression Ug, = + f/J as the full source term. We would, of course, prefer to 
incorporate some or all of into the definition of ([4.5|) in order to obtain a more compact 
expression for the source term. Unfortunately, as we shall soon argue, it would seem that 
this is not permissible, as doing so would destroy the desired equivalence with Einstein's 
theory. However, it is possible to re-express the combined source term in such a way that 
the residual torsion 2-form is not so manifestly apparent. Indeed, following the alternative 



route mentioned in the footnote before Eq. ( |4.8|) , one finds a combined source term which is 



almost entirely "torsion- free." We have chosen the route that we have, because it happens 
to lead to more aesthetically pleasing equations than does the alternative. That is to say, 
only exact exterior derivatives appear in the combined source term as we have written it. 
Following the other route, one finds a lot of triad covariant derivatives, -Db's, surfacing in the 
expressions. To our mind, the chosen route is more in spirit with the DF language intended 
to be showcased here. 



C. Wave-equation Formulation of General Relativity 

For the vanishing of Ta to produce a good wave-equation, the term Va. must be equal to 
a 2-form which involves fewer than second derivatives of H. As described by aacy, there 
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appears to be a number of ways to do this. For our purposes it suffices to enforce the simple 
harmonic slicing condition, 0| 

f = doN + N^H = 0. (4.14) 

(An obvious generahzation of this equation is obtained by assuming that / = f{t^x'^) is a 
prescribed well-behaved function.) One immediately sees that vanishes for a harmonic 
time-slicing. The vanishing of K, also removes the term involving triple derivatives of the 
lapse. Further, two terms in the expression ( [4. 11] ) for also vanish for harmonic time- 
slicing. 

As the new Einstein system of equations of motion, we now offer 

b^el = (e) 

□ tT, = -iV-^a (7f) 

doN = -N^H . (TV) 

The appropriate Cauchy data on the initial slice are e*, vTa, > 0, and -Do^Ta subject to the 
following constraints: 

= -47r[^^l (J) 
= -2Dn, (H) 
= 2VabVr^" - H7r\ + A /^a (H) 
= -2N-'bo7r^ + N-'el^ A D(A^a'^) - - H^e* + XcV^^e: + 2H7r, - , (tt) 

where the understanding is that these constraint equations are pulled back to the initial-data 
surface. Notice that we are requiring the full Einstein 3-form to vanish on the initial data 
surface. Both the shift vector and the skew-symmetric rotation matrix 0ab play a passive 
role in our formalism, as they have been absorbed into the operator Dq. 



V. VERIFICATION OF THE ORIGINAL EINSTEIN EQUATIONS 

Before showing that the new system of equations is equivalent to the original Einstein 
equations, we should demonstrate that it is indeed solvable. However, for sake of brevity, 
we only note here that the treatment of this issue found in the third reference of goes 
through essentially unaltered for the case at hand. Therefore, let us turn to the proof of 
equivalence. 

Claim : The vanishing of the Einstein 3-form and torsion is equivalent to the 2-form equation 

Ta = 0, (13) 

when subject to the initial-value conditions 



"Gx|s = 


(5.2a) 


"Gals = 


(Hd) 


Ga±|E = 


(B) 


T|e = 0. 





Here |s means "pull-back to the initial-data surface." 
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A. Method of Proof and Two Lemmas 



To prove the claim, it proves most convenient to switch to and proceed in index notation. 
Our proof follows the arguments of AACY and of York ||^ quite closely. We shall have need 
of the spacetime covariant derivative which "sees" quasi-coordinate indices, and we use 
to denote this covariant derivative. Consider the action of Sq and Vj on S scalars (such 
as Goo)) ^ 1-forms (such as Goi), and E 2- index tensors (such as Gij) which arise from the 
various components of a 2-index spacetime tensor (such as Ga/3) ■ The proof is facilitated by 
relating the action of do and Vj on these S objects to the action of Vq and Vj on the same 
objects. Such relations follow from the definition of these operators and the values of the 
connection coefficients associated with the quasi-coordinate frame (|1.1|), those being given 
by 

= Ig'^^iexlg.,] + eM - e^g^,]) + Ua'^^'g^uC^x" + g'^^gxuC,,'' - Cx,n ■ (5-3) 

Here the C^^^^ = [e^,ey]'^ are the structure functions associated with the frame (|1.1|) . The 
frame choice ( p. . 1|) ensures that the T^jk = T^jk are the ordinary Christoffel symbols associ- 
ated with the E coordinate frame. Throughout the proof, we shall employ the short-hand 
(l.h.E.T.)„ to mean "additive remainder terms linear and homogeneous in both the Einstein 
tensor and the torsion tensor and their derivatives of order < n." When remainder terms 
are built solely with the Einstein tensor or solely with the torsion tensor, we shall employ 
the short-hands (l.h.E.)^ and (l.h.T.)^, whose meanings should be clear. 

We begin by defining a S tensor Tij via Ta = eJe^'^Tije^. Provided that the triad is 
nowhere degenerate, the vanishing of Tij <^==^ the vanishing of Ta. Direct calculation now 
yields 

Tij = -doGij + WiGoj + V.Goi - h.VGoi - a, + (l.h.E.)o . (5.4) 

Hence, in ( |5.4| ) no derivatives of the Einstein tensor appear other than those explicitly shown. 
The S tensor Cij is defined via ej ^"^CijC^ = D{NTg), and for it we find the following explicit 
expression [cf. the footnote before 



= iV' {n^ [ViiN-%i) + V,iN-Xi) + ViiN-'T,^)] } . (5.5) 

The i ^ j symmetry present in the first two terms proves to be quite important for our 
construction. In fact, this symmetry (in tandem with the way the particular powers of the 
lapse enter the expression), ensures that the antisymmetric piece of N~^Cij is the Laplacian 
on Tij, up to a trivial (l.h.T.)g term. Compare the definition ( |5.4| ) of Tj^ with aacy's 
definition of Qij (or, even better, with their ^ij)- In this comparison, be sure to note that 
both doGij and Cij are not manifestly symmetric in our formalism. 

With the aforementioned relations between the actions of do and Vj and of Vo and Vj, 
we write 

T,, = -VoG,, + V.Go, + V.Go, - h,,V'Gok - C,, + (l.h.E.)^ ; (5.6) 
and, moreover, we have the following two lemmas for Cij. 
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Lemma 1: The expression for Cij may be written as 

= lN[V'V,Tji + VV.Ta + W^T,,-] + A,, + (l.h.E.T.)o , (5.7) 

where Ajj is a symmetric remainder term of (l.h.T.)^ order. 

To prove the lemma simply expand the expression (|5.5| ) for Cij. Then the proof reduces 
to verifying that VfcVzTjj = VfcViTjj + (l.h.E.T.)^. Now, viewed as a spacetime tensor, the 
torsion T^^ is tangential to the S slices (Tq, = Too = 0). It follows that VkTij = VkTif, and, 
hence, the only troublesome terms which arise in obtaining the verification above are those 
of the form VoT^. But such terms may be eliminated in favor of (l.h.E.T.)^ ones via ( p.9|) 
[technically, via ( p.9|) as we write it later in (|5.13| c)]. 



Lemma 2: The tensor Cij obeys the property V*Cjj = (l.h.T.)-^. 

Viewed as a spacetime tensor, C^js is by definition purely spatial (Coi = Qo = Coo = 0); 
it follows that V'^Cij = V^Cij. With this fact, one can prove the lemma directly by forming 
V^Cij and making appeals to the Ricci identity. Alternatively, with the definition of Cij one 
may show that 

ea^ V^a, = *D\NT^) . (5.8) 

Therefore, since Ta is (l.h.T.)^ by construction, the Cartan identity for the creation of 
curvature (symbolically = RA) yields the result (provided a non-degenerate triad). 



B. Zero-order Wave Equation for the Torsion 

Before turning to the main part of the proof, let us derive a certain zero-order wave 
equation for the residual torsion tensor Tjj which is implicitly defined by the vanishing of 
Tij. This derivation will also further elucidate the precise role played in our formalism by the 
1-form Ta appearing in the definition ( [4.5|) . The form of (|5.7| ) ensures that the antisymmetric 
piece of Cij obeys 

q,,] = iiVV'VzT,, + (l.h.E.T.)o , (5.9) 

Note that the extra terms in this equation are linear and homogeneous in the Einstein and 
torsion tensors with no spatial derivatives of them appearing. In particular, the tensor Xij 
has been killed by the antisymmetrization. 

Now consider the index version of (|2.9|), which may be expressed as 



VoT,, - 2iVG'[,,] = (l.h.T.)o . (5.10) 

With Eq. ( p.9|) in mind, we apply V° [the same as applying A^~^Vo for our frame (|1.1|) ] to 
equation ( ^.10|) and expand. Along the way, we make further substitutions with (|5.10|) to 
replace Vo(l.h.T.)Q terms with (l.h.E.T.)Q ones. We then reach 

V°VoT,, - 2Ar-iVoG[,,] = (l.h.E.T.)o . (5.11) 

Next, substituting ( ^.6|) into ( |5.11|) and appealing to (|5.9|) , we obtain 
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□T,, = -2N-^Tiij] + (l.h.E.T.)o , (5.12) 

with □ denoting the standard covariant D' Alembertian V^V^. Eq. ( ^.121 ) shows that the 
vanishing of Tij yields a zero-order wave equation for the residual torsion 2-form. This 
zero-order character proves to be quintessential in establishing equivalence between our DF 
system of equations given above and the vacuum Einstein equations. In fact, we have rigged 
the definition [cf. Eqs. ( |C5|) and (|C^)] of Ta in ( [4.5|) to ensure this zero-order character for 
( p.l2| ). Note that there is a projection into S implicit in ( p. 12 ). That is to say, while Tq^ = 



for the torsion viewed as a spacetime tensor T^^, it is not the case that DToj, the tangential- 
normal piece of DTq^, need be zero (of course, OTqq does indeed vanish by antisymmetry). 
Eq. (|5.12|) has nothing to say about DToj; however, for our purposes this observation is 
inconsequential. 



C. Twice-contracted Bianchi Identities and Equivalence Proof 



We shall require the so-called twice-contracted Bianchi identities V^G^^ = (index 
ordering is important!), which we expand out. 



= V^Goo + V^G.o 
= V^Go. + V^G,,. 



(5.13a) 
dnib) 



As is well-known, these identities are related to the preservation in time of the Hamiltonian 
and momentum constraints. Likewise, our identity ( |5.1CI| ), which we conveniently write now 
as 



= VoT,, + (l.h.E.T.)o 



is derived from the Bianchi identity for the torsion and is related to the preservation in time 
of the rotation constraint. 

As the first step in the proof, form V'Tjj, use the Ricci identity to commute covariant 
derivatives, insert ( ^.13| b), and appeal to the second lemma, in order to obtain 



□Go, = (l.h.E.T.)i. 



(5.14) 



As the second step, apply □ to the identity ( |5.13| a). Repeated use of the Ricci identity on 
the resulting expression yields 



V°nGoo + V^nG,o = (l.h.E.),. 



(5.15) 



At this point, we insert the result Gjo = Goj + TjkV^N [cf. Eq. ( 2.14 )]. Then, appeals to 
the earlier findings ( |5.12|) and ( p.l4|) show that 



VonGoo = (l.h.E.T)2. 



(5.16) 



Finally, as the proof's third step, apply □ to the definition of N^^Tij (which vanishes) as 
given in (|5.6| ), expand out the resulting expression, repeatedly use the Ricci identity, and 
make several appeals to ( [5. 14] ), thereby arriving at 
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N-^VoOG^j + a{N-^a,) = (l.h.E.T.)2 . (5.17) 

The lapse factor has been included to remove the lapse term in Cjj [cf. the form of this 
tensor given in the first lemma, Eq. ( |5.7|) ]. Note that n[N~^Cij) is (l.h.T.)^; however, that 
it is also (l.h.E.T.)2 may be shown by the following argument which makes use of ( |5.7| ). For 
the leading terms in N~^Cij, those quadratic in the covariant derivatives, we can move the 
action of the □ onto the torsion tensors in the expression. In so moving the □, we must 
again repeatedly use the Ricci identity, a process which generates only byproduct (l.h.T.)2 
terms. For n[N~^Xij), « priori (l.h.T.)3, we note that, again up to (l.h.T.)2 terms, the □ 
can be moved onto the torsion tensors in the expression. Therefore, invoking ( |5.12| ), we get 

VoOG,, = (l.h.E.T.)2. (5.18) 

Notice that the application of □ has mapped (l.h.T)2 terms into (l.h.E.T)2 ones. The ability 
to take this crucial step rests squarely on the ^ero-order character of the wave equation 

With ( |5.14| ), ( |5.1(j| ), ( |5.18| ), we obtain the following system of equations: 



doaGoj = (l.h.E.T.)2 (5.19a) 
doDGoo = (l.h.E.T.)2 

doOG,, = (l.h.E.T.)2 ( CT c) 

4t,, = (l.h.E.T.)o. (lO^ d) 

This system is (up to homogeneous Einstein-tensor and torsion-tensor terms) essentially 
the system considered in the works of AACY, although here we see the novel appearance of 
the torsion tensor. It is a strictly hyperbolic system. All (l.h.E.T.)2 and (l.h.E.T.)o terms 
on the right-hand side may be shown to be zero on the initial Cauchy surface by (i) the 
(index versions of the) boundary conditions (|5.2|), (ii) the Bianchi identities ( p. 13] ), (iii) the 
vanishing of Tjj, and (iv) derivatives of the previous three items. We infer that Tij and all 
pieces of the Einstein tensor save Gjo vanish in the future domain of dependence D~^(T,). 
We then get the vanishing of Gjo from the vanishing of Goj and Tij by the identity ( p. 141 ). 
Equivalence has thus been established. 



VI. CONCLUSION 

We conclude by giving some basic instructions (with cautionary remarks) for rewriting 
our DF results in the tensor-index language. At the outset, we must admit that we find 
passing back and forth between the DF and tensor-index versions of wave-equation general 
relativity to be not as clean and direct as we had originally hoped would be the case. Passage 
to the tensor-index version of the formalism essentially amounts to (i) defining Tj^ as before 
in Eq. ( p.4[ ) and then (ii) simply setting the torsion to zero by hand in all expressions [hence 
Cij = in (|5.4| )], thereby obtaining a symmetric Tij. Now, it is tempting to simply use Tjj 
as the jumping-off point for deriving an tensor-index wave equation for TTjj, working with Tij 
in the same fashion that AACY work with their Qij. However, as we shall now argue, though 
this is essentially the correct procedure, one must be careful to extract superfluous Einstein 
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tensor terms which are indeed present in Tij. In ( p. 41) these superfluous terms have simply 
been swept into the (l.h.E.)Q short-hand, but let us now examine the origin and nature of 
these terms in more detail. 

Consider again the source term ([4 .111) written as follows: = ej e^'^Uf^e^. Now, in fact. 



f/J still has some Einstein tensor "buried in it," i. e. one flnds that f/JI- contains contracted 
covariant derivatives of vTjj, e. g. Viir^, coming from the sixth term in the expression for U^. 
But 2V;7r'j = G±i, and, hence, in the 3+1 decomposition of Tij superfluous terms involving 
the momentum constraint G±i arise. Therefore, one should locate and remove these super- 
fluous momentum-constraint terms in the passage to the tensor-index formalism. However, 
besides such superfluous momentum-constraint terms, another superfluous Einstein-tensor 
term also arises in the passage to the tensor-index version of the formalism. Indeed, consider 
the leading term in Ta, namely, 

^oGa± = -eJe''\doG,,)el - G,,bo{e^' e'^'el) . (6.1) 

The expansion ( |S.1| ) shows that our E tensor Tjj has an undifferentiated Gij term "buried" 
in it; hence, in passing to the tensor- index language we be should be careful to locate and 
remove this term. 

One might ask, why not simply work with a modifled deflnition of Ta (call it T^) such 
that passage to the tensor-index formalism is rather more direct, i. e. no Einstein tensors 
surface in the 3 + 1 decomposition of T[- (obtained from T'^ = Ca-'e^^T-^e^). However, 
the differential-form wave equation for vTa obtained from the corresponding 2-form equation 
(T') is quite a bit messier than the one we have considered, and, if fact, this wave equation 
contains superfluous Einstein 3-form terms in it! Therefore, if the intent is to remain within 
the framework of DF, then we flnd it preferable to work with ( |5.1| ) and its 3+1 decomposition 
as given. However, this means that one must not be too naive when rewriting the results 
in the tensor-index language. One must be careful to extract superfluous Einstein-tensor 
terms which will arise in the tensor equations directly obtained from the corresponding 
DF equations. The non-triviality associated with going back and forth between the two 
languages can be traced to whether it is the Einstein 3-form or the Einstein tensor which is 
held as the fundamental object (i. e. which is derived from which). 
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APPENDIX A: GAUSS-CODAZZI-MAINARDI EQUATIONS WITH TORSION 



The GCM equations (see for example [|l],|T2[ ) are integrability criteria which relate the A4 
Riemann tensor 3?"/3^,y to the S Riemann tensor R^jki, the E extrinsic curvature tensor Kij, 
and other geometric objects associated with the S foliation. The form of these equations 
is complicated quite a bit by the presence of spacetime torsion. We shall present equations 
valid only for the limited type of spacetime torsion considered in this work. 

To derive the needed equations, one may proceed in a number of ways. Perhaps the best 
method is to apply the form decompositions ( |1.3| ) and (|1.4|) to the second Cartan structure 
equation for 3?"^^, the curvature 2- form. This straightforward method leads to the following 
expressions for the time-gauge components of the Ai Riemann tensor: 

3^abcd = -Rabcd + -^ac-^bd " -^ad-^'bc (Ala) 

3?±abc = D^K^ - D,K,^ (^) 
3f^ab±c = N~'^ec^douJs,hj - a^^Khc + a^K^c - N''^ (ec[rabo] - Tdbot^'^ ac - Tado'^'^ be) (^c) 
3^a±b± = N-^e^WoK^j + D^a, + a^a^ . (^) 



Let us express the identities ( [Al|) in terms of E coordinate indices. Since all of the terms 
in ( [Al^ ,b) are tensorial, conversion of these equations into the required forms is immediate. 
Moreover, with the evolution rule (|1.12|) for the triad, conversion of the final identity (|Al|d) 
into E coordinate indices is also not difficult. However, finding a suitable expression for 
[T|c) is more difficult. One can use the (PT^)^ piece (|2.10 ) of the Bianchi identity (|2.5| ) for 



the torsion 2-form, in order to express KabXc as ^±ca.h plus terms built from the torsion tensor 
Tab and the lapse function N (recall that without torsion the Riemann tensor is symmetric 
under such pair exchange). We obtain the following list: 

^ijki = Rijki + KikKji - KiiKjk (A2a) 

^±ijk = VjK,k - VkK.,j (Pb) 

= V,K,k - VkK,, + \N-^V^{NTjk) + \NVj{N-^Tk^) - \NVk{N-%^) (^c) 

3?i±j± = N-'doKij + K\Kij + Vjfli + aiaj , (Pd) 

where we remind the reader that Tij = 2K[ij], and Vj stands for the ordinary E covariant 
derivative. Clearly, equations (|A2| b,c) exhibit the loss, due to torsion, of the pair-exchange 
symmetry enjoyed by Riemann curvature in the absence of torsion. However, note that 



(|A^d) is also not symmetric under i ^ j exchange. However, all of the terms on the right- 
hand side of (|A^d), save the very first, are symmetric under i ^ j exchange. In particular, 
note that the term VjOj = VjVjlogA^ is indeed symmetric, as we have set the E torsion 
tensor T^jk manifestly equal to zero. The behavior of = ^±i±j under the pair exchange 

±z ^ ±j is also evident from (|2.7|) . 



APPENDIX B: COMMUTATOR OF Dq AND D 

To evaluate the commutator [Dq, D]ip^, where ip^ = ^kp^ is an arbitrary triad-valued 
spatial form, it suffices to consider the action of the squared spacetime exterior covariant 
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derivative on an arbitrary tetrad- valued spatial form ip"^ = i'?/'^. Applying the form 
decompositions (|1 .31 , 11 .41) to the expression Vip"^, one finds that 

\ViIj^) = diP^ -K^A-^^ (Bla) 

{Vij^)± = N-'do^fj^ + a^^fj"" (il|b) 

\V^^) = -K^A^^ (|B3c) 

(p^'^)^ = N-'^DoiP'' + a^'ip^ . (^d) 

Now, a straightforward expansion yields that 

p2^^ = D^Vip^) -K^ A ^Vip^) 

A [N-^[Do\Vtlj'') - DNiVij"")^] + a^'^V^^) + K"" A . (B2) 

However, by the second Cartan structure equation we know that the above expression is also 
p2^a — s)fja ^ y\ where as before 3fJ"^ b is the matrix- valued curvature two-form of Cartan, 
which in index notation is represented by the mixed-components '^^b^u of the spacetime 
Riemann tensor. Therefore, taking the i\_ inner derivative of (|B2|) , inserting the relations 
and collecting like terms in ip-^ and we determine the following identity: 

(W" x± + N-^DoK"" + Da"" + a^a) Aip^ = 

N-\DoD - DDoW - {^^ h± + a'^K^ - a^K'') A tp'^ . (B3) 

Since each side of the above equation is independent of the other, we get two relations. The 
first yields the "dynamical" GCM equation ([Al| d). The other equation we get from ( |B3| ) is 
the desired one (|1.14|) . 

In the expression (|1.14|) the curvature term can be eliminated in favor of extrinsic- 
curvature and acceleration (ob = -DblogA^) terms. Indeed, using the GCM equation (|A^), 
one may derive the following identity: 

3?ab±d + 2a[,irb]d = iV-i[D,(iVir(db)) + /^d(iVi^(ab)) - /^b(iVi^(da)) - ^d(iVi^ba)] • (B4) 



Remarkably, on the way to this result all of the torsion terms stemming from ([A^c) conspire 
in such a way that only the symmetric piece of the extrinsic curvature tensor appears in 
all but the last term above. If desired, this identity may be inserted into the result ( |1.14D 
for the commutator. As an application of these considerations, consider a triad-valued 0- 
form /b = eb*/i obtained by "soldering" the free index of a S 1-form /j. With (p^) and 
the evolution rules ( |1.12| ) in hand, a careful examination of the expression j ig. ([l)o, -D]/b) 
(note the "hook" with di) yields the tensor- index form of the commutator. Namely, 

[4, Vi]/, = [Vi{NK(,k)) + V,{NK^,k)) - Vk{NK^,,))]f . (B5) 

This commutator, expanded out and written in terms of the acceleration 1-form Oj = 
Vi log A^, is featured prominently in the works of AACY. 
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APPENDIX C: TIME DERIVATIVE OF THE TRIAD CURVATURE 



Recall that in index notation we have the following definition of the triad curvature 2- 
form: R^ij = —^€ahcR^^ij, where the R^^ ij are the mixed triad-coordinate components of 
the S Riemann tensor. We begin by using the GCM equation (|Al|c) to derive the result 



DoR. = -ieabc/^(iV3fJ''' ± + 2Na!'K') . (CI) 

This identity is essentially the {VW^ b)± piece of the full (or "un-contracted" ) Bianchi identity 
Vm^B = for the spacetime curvature 2-form. Indeed, another method to derive ( |(J1D 
involves a 3 + 1 form decomposition of the full Bianchi identity. 

We now derive several auxiliary results which will be used with ( PT| ) to get the appro- 
priate expression for DqRs.- First, using the exphcit expression ( p.6|) for vTa along with the 
GCM equation (|Al| b), we obtain 

|(3?xa"^)e^bc = -^'D\^ - {D'^H)el^ " (^'T^)e^, . (C2) 

Second, performing some permutation-symbol gymnastics and appealing to the identity 
(^) we find 

|eabc3fi"\ - |(3fJ^\a)eL = 'C'^^'X^, + CabcT^a'^e'^ . (C3) 



Third, with the definition ( ^.81 ) we straightforwardly write 

eabcfl'i^^ = -2eabca' V + Ha'^el^ - e,^,T\a'e'' . (C4) 

Notice the c ^ d antisymmetry present in the last terms of ( |C3| ) and ( |C^ ) [we have not 
employed the e*^^ notation in (|C^ ) in order to highlight this antisymmetry]. Finally, we 
define the triad-valued 1-form Ta with the seemingly strange equation 

- i(3fJ±a"^ - 3fJ'^±a)e^, = Ta + (z^^T^ - a^^T^a + iflaT" j el, . (C5) 



We have thus defined Ta, because we demand that the tensor Cij ( [5.5| ) stemming from Ta 
satisfies the property (|5.9| ). This property would seem essential for the proof of equivalence 
given in § V to go through. Although it may yet be possible to modify the definition of 
Ta in an advantageous way, the definition that we adopt here provides us with a consistent 
formalism. Using the GCM splitting equations (|A^,c), we find^ 

Ta = i5abiVrc^(iV-'T)]e" - 6^kN*D\N-'e'' A ^T) . (C6) 

as the abstract expression for this 1-form. 



^Alternatively, we can write 



-k 

ebc 



This results allows one to easily obtain the explicit expression for Q,- given in Eq. (^?5| 



24 



Now add together the equations ( |C2| ) , ([031), (|C^), and (^) and then perform a short 
calculation (taking advantage of the aforementioned c ^ d antisymmetry) in order to reach 
the following result: 



a 



-2eabca'^ + (a^i/ - I^''i^)e:b + a^'T . (C7) 
Finally, plugging (|C^ ) into (|Clj ) we arrive at 



boR^ = -iV'D^iiGa - e*b A D*(iV"G^) - D{NT^) - 2NAtt^ + 2Na A *D*tt^ 

+2eabc^(A^a^ V^) + e^^ A D{NHa^) - e^^ A D{ND^H) - D{Na^*T) . (C8) 

To get this identity, we have introduced the triad-covariant Laplacian A = {*D*D — D*D*) on 



triad- valued 2- forms at the price of an appeal to ( |3.14| a). Also, in this expression we consider 
T as a short-hand for 2e^ A [cf- Eq. ( |3.10|) ]. We remark that other useful expressions 
for the final form of -Do-Ra are certainly possible. However, although we would, of course, 
prefer to incorporate the last torsional term appearing in (|C8|) into the definition of Ta (in 
order to obtain a more compact expression for the final result), this would seem not to be 
permissible. Indeed, were we to alter the definition of Ta in this manner, we would lose 
the crucial condition (p.9|). More precisely, we would have to replace the (l.h.E.T.)o in that 
equation with (l.h.E.T.)i. This in turn would seem to spoil the proof of equivalence between 
our DF system of equations and the original Einstein equations. 
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